Finite element approximation of one-sided Stefan problems with anisotropic, approximately crystalline, Gibbs-Thomson law
نویسندگان
چکیده
We present a finite element approximation for the one-sided Stefan problem and the one-sided Mullins–Sekerka problem, respectively. The problems feature a fully anisotropic Gibbs–Thomson law, as well as kinetic undercooling. Our approximation, which couples a parametric approximation of the moving boundary with a finite element approximation of the bulk quantities, can be shown to satisfy a stability bound, and it enjoys very good mesh properties which means that no mesh smoothing is necessary in practice. In our numerical computations we concentrate on the simulation of snow crystal growth. On choosing realistic physical parameters, we are able to produce several distinctive types of snow crystal morphologies. In particular, facet breaking in approximately crystalline evolutions can be observed.
منابع مشابه
Approximation of Crystalline Dendrite Growth in Two Space Dimensions
Abstract. The phase transition between solid and liquid in an undercooled liquid leads to dendritic growth of the solid phase. The problem is modelled by the Stefan problem with a modified Gibbs-Thomson law, which includes the anisotropic mean curvature corresponding to a surface energy that depends on the direction of the interface normal. A finite element method for discretization of the Stef...
متن کاملApproximation of Crystalline Dendritegrowth in Two Space
The phase transition between solid and liquid in an undercooled liquid leads to dendritic growth of the solid phase. The problem is modelled by the Stefan problem with a modiied Gibbs-Thomson law, which includes the anisotropic mean curvature corresponding to a surface energy that depends on the direction of the interface normal. A nite element method for discretization of the Stefan problem is...
متن کاملStable phase field approximations of anisotropic solidification
We introduce unconditionally stable finite element approximations for a phase field model for solidification, which take highly anisotropic surface energy and kinetic effects into account. We hence approximate Stefan problems with anisotropic Gibbs–Thomson law with kinetic undercooling, and quasi-static variants thereof. The phase field model is given by θwt + λ %(φ)φt = ∇ . (b(φ)∇w) , cΨ a α %...
متن کاملOn stable parametric finite element methods for the Stefan problem and the Mullins-Sekerka problem with applications to dendritic growth
We introduce a parametric finite element approximation for the Stefan problem with the Gibbs–Thomson law and kinetic undercooling, which mimics the underlying energy structure of the problem. The proposed method is also applicable to certain quasi-stationary variants, such as the Mullins–Sekerka problem. In addition, fully anisotropic energies are easily handled. The approximation has good mesh...
متن کاملExistence of weak solutions for the Stefan problem with anisotropic Gibbs-Thomson law
The Stefan problem with Gibbs-Thomson law describes solidification phenomena for pure substances. In applications the surface energy is anisotropic leading to an anisotropic Gibbs-Thomson law. We show the existence of weak solutions to the Stefan problem with anisotropic Gibbs-Thomson law using an implicit time discretization, and variational methods in an anisotropic BV setting. Our main resul...
متن کامل